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Abstract 

We show that in the Two Higgs Doublet Model (THDM) Type I the partial decay width 
T(h — > 77) can be considerably larger than in Type II, due to light charged scalars with 
M H ± 100 GeV, which are not yet excluded in Type I. A possible enhancement of the width 
compared to the SM is analyzed for different Higgs potentials, subject to constraints from tree- 
level unitarity, vacuum stability and electroweak precision tests. 



1 Introduction 

While most aspects of the Standard Model (SM) have been confirmed by experiment, the elec- 
troweak symmetry breaking mechanism has yet to be established. The SM realizes symmetry 
breaking with a single Higgs doublet, giving rise to only one physical neutral scalar. 

In the THDM one increases the scalar particle content by adding a second doublet, resulting 
in 5 Higgs bosons (h, H, A, H^). The most general version of this model is however problematic, 
because it induces Flavour Changing Neutral Currents (FCNC) in the Yukawa sector. One usually 
introduces discrete symmetries to solve this problem, leading to the Type I and Type II versions 
of the model [l|. Having the Yukawa structure of the Minimal Supersymmetric Standard Model 
(MSSM), the Type II model has been studied extensively in the literature. Despite its interesting 
phenomenology, the Type I model has received less attention in the literature and has mostly been 
studied in the fermiophobic limit @-Q|- 

In the Type II model one can derive a strong exclusion limit for the mass of the charged scalar 
M H ± > 300 GeV independently of tan (3 (ratio of scalar vacuum expectation values) because of 
indirect constraints from B — > X s j [H]. Due to the different Yukawa interactions, the situation for 
M H ± is different in the Type I model. The strongest lower bound arises from searches at LEP, 
still allowing for a light charged Higgs with M H ± > 90 GeV [6]. As most of the early searches at 
the Large Hadron Collider (LHC) for Higgs bosons will focus on a SM-like light neutral Higgs, it is 
important to study how much deviation from the SM one can expect for the THDM. One suitable 
decay channel to search for such deviations is H — > 77, for which a(pp — > H)B{H — > 77) may be 
measured with a relative error of 10 to 15% for 100 fb _1 integrated luminosity for a SM Higgs H 
with 110 GeV < Mjj < 150 GeV 0- 

The decay h — > 77 was studied previously in Type II in Refs. [1, @], which showed that an 
enhancement of the partial decay width T(/i — > 77) by around 25% may occur. 
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For Type I, this decay was studied in the fermiophobic limit in Refs. 0, In this limit, 
B(h — > 77) might get considerably enhanced because h — > bb is loop-suppressed, resulting in a 
smaller total width of the Higgs boson h compared to the width of the SM Higgs. Such fermiophobic 
Higgs bosons might hence lead to improved detection scenarios via h — > 77 at the LHC 0]. 

On the other hand, cancellation effects may lead to V(h —> 77) ~ 0, which would render Higgs 



detection via this decay impossible [U] 



In our study, we are interested in enhancement effects of T(h — > 77) itself. While —> 77) 
can be measured very precisely at the International Linear Collider (ILC), enhancement effects of 
it may also be visible early on at the LHC if the Higgs production in the THDM is not suppressed 
compared to the SM. Recently, such effects were studied in the context of Higgs production at 
a photon collider In our study, we analyze the effect of light charged Higgs bosons with 

M H ± « 100 GeV for T(h -»• 77) in the THDM Type I in the range 110 GeV < M h < 150 GeV, 
constraining the parameter space with constraints from vacuum stability, unitarity, and electroweak 
precision tests. The constraints considered are more restrictive than the ones of Ref. and 
therefore the possible enhancement we find is considerably smaller but yet still enough to be seen 
at the LHC. 



2 Background 

After introducing the Z2 symmetry $1 — > $1, $2 — > —&2 for the two Higgs doublets $1, $2 arid 
allowing for a soft breaking term, we get the potential 



V = m\ x <&[<$>i + m2 2 3>2*2 — x$2 + h.c 
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As we will not consider CP violation, all parameters are assumed to be real. We define the special 
cases 

• Va with m\ 2 = 

• Vb with A5 = 0. 

The naming convention follows Ref. [13|. After minimization of the potential , one introduces the 
vevs of the Higgs doublets $1, <&2'- 
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As we only consider minima which do not violate CP, we assume that both v\ and V2 are positive 
real parameters. Only the potential Va allows for either v\ or V2 to be exactly zero. In this case, 
the Higgs doublet with the non-zero vacuum expectation value must couple to the SM particles. 
This version is called the Inert Doublet Model (IDM) in the literature, and we will take v\ = in 
this model 14 ]. As the Z2 symmetry is unbroken in this model, the lightest particle of <I>i cannot 
decay and will contribute to Dark Matter. 
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In all versions of the model, v = \Jv\ + v\ ~ 246 GeV to yield the correct masses of Mz and 
Myy. As all parameters are real we can diagonalize the mass matrix of the two CP-even Higgs 
bosons with a single parameter cos a (see [l] for more details). 

For the potential V we can express all parameters but m\ 2 with the four scalar masses (M^, 
Mjj, Ma, M h ±), tan/3 and cos a. As the potentials Va and Vb have one parameter less, they can 
be completely expressed with the Higgs masses, tan f3 and cos a. In the IDM there are no mixing 
angles, and 2 parameters (Ai and m 2 ^) will remain unexpressed in this model. 

The Yukawa interaction in the Type I model are defined so that only <3?2 couples to the fermions. 
In our convention, the coupling hbb oc cos a/ sin /3, where h denotes the lightest CP-even Higgs boson 
(or the SM-like Higgs boson in the case of the IDM) . 



3 The constraints 
3.1 Vacuum stability 

If we want our vacuum to be not only a local, but also the global minimum, one has to impose the 
following conditions on the parameters of the potential in (pQ) [15]: 

Ai, A 2 > and - (AiA 2 ) 1/2 < A 3 

-(AiAs) 1 / 2 < A 3 + A 4 -|A 5 | • (3) 

These conditions yield important constraints on the Higgs masses and mixing angles. 



3.2 Unitarity 

Unitarity constraints in the THDM for the potential V were calculated in Ref. [16]. In this analysis 
various processes were used to constrain quartic couplings at the tree level. These lead to a set 
of unitarity constrained parameters a±,b±,c±,d±, f±,e± : 2, f 1,2, Pi, whose absolute values must be 
smaller than 8ir. These depend on the parameters of the potential, which in the case of V are the 
4 Higgs masses, tan/3, cos a and m 2 J^]. For the potential Va one simply sets m\ 2 = 0, and for Vb 
m 2 2 = M\ sin/3 cos j3. In the IDM the above constraints depend on the 4 Higgs masses, Ai and 



3.3 Ar 

Another powerful constraint comes from the well measured constant Gf = 1.16637(1) x 10~ 5 GeV~ 2 , 
which is defined via the muon decay fi — > e Vnp e (p/) in the effective Fermi theory. Calculating these 
decays in the THDM one can relate Gf to the self energies of the vector bosons [l8( ] : 

Gf a f , 

-^ = — 2 2- l + Ar . 4 

Here s\y is the sine of the weak mixing angle syy = sin 6\y with 



sin 2 ^ = l-^. (5) 
m z 



1 For a translation of the parameters used in Ref. [lg] to the ones employed in Eq. JT]) see Ref. [17 
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The quantity Ar is then separated into the finite quantities 



Ar = Aa - ^fAp + Ar rem , (6) 



s w 

^_E«-(0)_ 2£!I E^0) i AQS _ flen .« (m|) + n .« (0) . (7) 
m| mfo cw rn z 

4S xz (0)-fieS zz (m 2 7 ) a / 7 - 4sL , 2 \ 
+ _Hi U L^Z + 6 + — Jlogc^ . (8) 

*w "4 i7rs w V y 

Here Y, zz , E^" 2 denote the transversal parts of the self-energies of the and Z boson. II" 4 ^ 

is defined as 

u AA {k2) s , (9) 

where S" 4 " 4 is the transversal self-energy of the photon. Ap and Ar rem at the 1-loop level were 
calculated in FeynArts/FormCalc [lil. |20| and then evaluated with LoopTools 20, 21]. We take 
Aa(m|) = 0.0594(5), where most of the uncertainties come from the hadronic contributions 22j] 



Note that due to the Appelquist-Carazzone theorem all heavy particles like the top and the Higgs 
bosons decouple from Aa [23(. On the contrary, Ap and Ar rem are sensitive to the masses of the 
Higgs bosons, the W and Z bosons and the top quark. 

As Ap depends quadratically on the top mass, while Ar rem only depends logarithmically on it, 
we include 3-loop QCD corrections only for the Ap parameter, which were calculated in Rcf. [24] . 
If one takes into account 2-loop contributions in the SM, the value of Ar shifts up by ~ 0.005 [251 ] . 
We take this as a rough estimate for the error arising from neglecting the 2-loop contributions in 
the THDM. 

On the other hand, one can simply insert the precisely measured values mz = 91.1876(21) GeV, 
a(0) = 1/137.03599976(50) and m w = 80.398(25) GeV into flD, which then yields an experimental 
prediction for Ar: 

Ar exp = 0.0343 ± 0.0020 . (10) 

Combining the uncertainty from experiment at 2 a with the theoretical error due to neglecting the 
2-loop contributions, we can exclude Ar < 0.0253 and Ar > 0.0433. At the one-loop level Ar is 
only sensitive to 5 unknown parameters, sin(/3 — a), M H ±, Ma, M^, Mjj. For the top quark we 
use the value m t = 171.2(2.3) 0- 

3.4 Z decays 

Yet another observable which has been determined with great precision at LEP and SLAC is the 
hadronic branching ratio of Z to bb pairs: 

r(z -^bb) 

6= r(Z^hadrons)' ( ' 
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A recent data analysis in Ref. 



27] showed that 



R b = 0.2163 ± 0.0007 , (12) 
which is in good agreement with the SM, which yields Rf M = 0.2158 [i^]. The ratio was calculated 



in the THDM in Refs. 29] and 30]. For Type I, one finds that the contribution of the charged 
Higgs bosons dominates the contribution of the neutral Higgs bosons, and therefore we can neglect 
the latter. One can hence derive a lower bound for tan /3 that depends on the charged Higgs mass. 



Using the appropriate formula of [30] one finds that for rat = 171.2 GeV 

tan < 2 excluded for M H ± « 100 GeV (13) 

at 2 a. Increasing the Higgs mass the bound drops (e.g. for M^± ~ 500 GeV one can exclude 
tan/3 < 1). One can find similar bounds by considering B — > X s j |3l| H. 



4 Analysis of h — > 77 



Having introduced the constraints, we will apply them to constrain possible enhancement effects 
in T(h — > 77). The partial decay width T(h — > 77) is loop induced, and can be easily calculated in 
FeynArts. The 1-loop analytical result is well known (see [12] and references therein): 



T(h -> 77) 



G F o?Ml 



128\/2V 



^2 9fQ 2 fN c A l/2 {Tf) + g w Ai(T W ) + g h A (r H i 
fe{t,b} 



where 



gw = sin(/3 — a) 



9f 



cos a 



sin/3 



9h 



raw 
2 — 9hH+H~ 

gm H ± 



(14) 



(15) 



with the contribution of the other light fermions neglected, gw and gf are the trilinear couplings 
of h to the W gauge bosons and to the fermions normalized to the ones of the SM. ghH+H- 1S 
the coupling which appears in the Lagrangian, C = g hH + H ~hH + H~ ■ ■ ■ . 7 
factors Ai are defined as 



t%, and the form 

4m. ' 



Mr) 

Al/2(r) 

Mr) 



-[r-f(r)]r-\ 

2[ r + ( T _i)/( T )] T -2 , 

-[2t 2 + 3t + 3(2t-1)/(t)]t- 2 



and 



fij) 



arcsm 
1 

~4 



log 



1 + vT 



1-VT 



ITT 



T < 1 

r > 1 



(16) 



(17) 



Note that gf is the same for b and t quarks in Type I. To get the result in the SM, one simply sets 
9f = 9w = 1 and g^ = 0. The result for Type II is almost the same, only g t = cos aj sin/3 and 



2 The exact value of the bound is less important for our study. What is important, is that we get a lower bound 
for tan j3 which is well above 1 for M H ± w 100 GeV. 
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Rjj — . r . (18) 



gb = sin a I cos j3. If the coupling hbb is not enhanced, the contribution of the b quarks is negligible 
and the result is approximately the same for both types. The main difference between the Type 
I and Type II models hence comes from the different constraints available, particularly the light 
charged Higgs bosons, which are not yet excluded in Type I. As we only consider Type I, the 
contribution of the b quarks (and the other fermions) will be neglected in the following discussion. 

In the parameter space region where 110 GeV < Mh < 150 GeV and Mjj± = 100 GeV, we can 
see that the contribution of the W bosons is by far dominant, being around — 10 < A\ < —8, while 
0-4 < A) < 0.5 and A 1/2 « lM Due to the fact that gf may only rise for small tan/3, which is 
restricted by the constraint, the most significant enhancement arises from g^. 

For a useful comparison we define the ratio i? 77 , the partial decay width normalized to its value 
in the SM: 

T(h -> 77) 
T(h -)• 77)5m' 

We then maximize this ratio subject to the constraints introduced in the previous sections. For the 
maximization procedure we assume 110 GeV < Mh < 150 GeV, and Mjj± = 100 GeV for the mass 
of the charged scalars. For the masses of the other Higgs bosons we assume Mh < Mh < 1 TeV 
(except in the IDM, where we may also have Mh > Mh), and 100 GeV < Ma < 1 TeV. The 
maximization was done in Mathematica 6, using its three different maximization algorithms. While 
these algorithms do not prove that a certain point is a global maximum, they easily find the regions 
where the width is enhanced. 

4.1 IDM 

As the Higgs boson h of the IDM behaves exactly like the SM Higgs when coupling to SM particles, 
we have gw = 9f = 1 (which corresponds to a = 0, /3 = 7r/2). For the trilinear coupling to the 
charged Higgs bosons we get ghH+H- = — A3U, which results in 

M 2 ± - ml. 

As Ai is negative and Aq positive, we must have m\ 1 3> M H± to get constructive interference. This 
means we must have A3 < and —A3 3> 1, which is constrained by the vacuum stability conditions 
in Eq. ([3]). If mf 1 i.e. A3 is large, Ai must also be large to be compatible with the constraints 
from vacuum stability (A2 cannot compensate for Ai, because we have A2 = M^/v 2 and Mh is fixed 
in our discussion). Furthermore, large values of Ai and A3 are constrained by unitarity, especially 
|o+| < 8tt. We hence find maximal enhancement to be around 1.6 < i? 77 < 1.8 in the region 
where 110 GeV < M h < 150 GeV for light charged scalars with M H ± = 100 GeV. As M 2 H± < v 2 , 
increasing M H ± will not decrease a_|_ significantly to allow for larger m 2 ^, and only results in a 
smaller gh and therefore an overall smaller enhancement. We plot a possible enhancement region 



in Fig. 1(a) As shown in Fig. 1(b) , further increasing Ai weakens the vacuum stability constraint, 



but increases the region forbidden by the unitarity constraints. 

If mfj is negative, the contributions interfere destructively. The contribution of the charged 
Higgs bosons must hence get more than twice the contribution of the W bosons to get enhanced. 
Using \a + \ < 8tt one can derive the rough bound mn < 600 GeV, which yields i? 77 < 0.9, and 
therefore no enhancement is possible in this case. 



3 For plots of these functions, see [32 
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100 150 200 250 300 350 100 150 200 250 300 350 



(a) i? 77 for Ai = 8.2. The allowed region increases (b) _R 77 for Ai = 8.355. We see that unitarity 

slightly for higher values of m^. constraints forbid larger values of Ai for high values 

of mn. 



Figure 1: Curves of constant value of it! 77 are plotted in the (M^, mu) plane for the IDM. The 
exclusion regions were calculated for M H ± = 100 GeV, Mjj = 90 GeV and Ma = 120 GeV. 
Ar ~ 0.038 in both regions displayed. 



4.2 V A 

For the potential Va (with both v\ , V2 ^ 0) , we have 

9h = 91+92 , 



(20) 



where 



91 = S/3- 



2Ml ± 



92 



Ml 



S2P M 2 H± 



(21) 



As gi oc g w , a significant enhancement is only possible for large M^. In the region we consider, 
—0.13 < gi < 0.4 and therefore it cannot significantly contribute to an enhancement of g^. 

g2 enhances the contribution of the charged Higgs bosons for large tan/3. To make the interfer- 
ence effects with the W bosons more obvious, one rewrites Eq. (|14p as 



T(h -> 77) 



G F a 2 Ml 



128\/2- 
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where 



C = spAi + 



A* 



Ml 



sp 2M 2 
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CfiA x + 



Ar 
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c p 2M 2 H± 



(22) 



(23) 
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with A 



1/2 



QtN c Ai/ 2 (T~t) and the term proportional to g\ neglected. As both C and D are real 



for the parameter region considered, maximizing over c a , s a simply yields 



max {T(h — > 77)} 



G F o?Ml 



128^/271-3 



C 2 + D< 



(24) 



Since Aq and A\ have opposite sign, the contributions of the charged Higgs and the W bosons 
interfere destructively in both C and D. 

Expanding C and D for large tan/3, one sees that C is of order 1, while in D the contribution 
of the W bosons is of order 1/tan /3, and the one of the charged Higgs bosons is of order tan /3. For 
moderate values of tan (5 (i.e. 3 < tan /3 < 6), C is by far dominant, and D is only around 5% of C. 

When raising tan/3, \D\ first drops to zero as the two contributions cancel, and then starts to 
grow again as the contribution of the charged Higgs bosons starts to dominate in D. One would 
hence expect an enhancement from the point on where D compensates for the suppression in C, 
which happens at tan /3 > 8. However, when employing the maximization procedure we find that 
such large values of tan/3 are strongly restricted by unitarity (especially a + ), and only moderate 
values of tan /3 are allowed. The partial decay width hence gets reduced compared to the SM, with 
typical maximal values -R 77 ~ 0.8 for 110 GeV < < 150 GeV as shown in Fig. [2j As we stay 
in the region where C dominates, raising Mh raises the negative contribution of the charged Higgs 
bosons in C and leads to a slightly stronger suppression compared to the SM. 



Unitarity violated 




Xanfi £ 2 excluded 



-0.5 



0.0 



0.5 



03 



Unitarity violated 




tan/? % 2 excluded 



-0.5 



0.0 



0.5 



(a) Ry-y is shown for M h = 110 GeV. Ar sa 0.037 
in the region displayed. 



(b) i? 77 is shown for Mh — 150 GeV. Increasing 
Mh decreases C and hence reduces _R 77 . Ar ~ 
0.038 in the displayed region. 



Figure 2: i? 77 analyzed for the V A potential, for M H ± = 100 GeV, M A = 100 GeV, M H = 150 GeV. 
We see that there is no enhancement, but a reduction compared to the SM. The displayed region 
is not forbidden by Ar or the vacuum stability conditions. 
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4.3 V and V B 



At last we will study the potential V (and Vb), where mf 2 is a f ree parameter. We now have 

9h = 91 +92 + 93, (25) 

where g\ and g<i are again the parameters defined in (|2ip . g% is not fixed via the masses, but is 
proportional to 



m 



.2 . 
12- 



.93 



(26) 



^3 has the opposite sign of #2- The interference behavior for V is again illustrated by Eq. (G 
where now in C and D we must replace 



Ml 
2M 2 H± 



M 2 h 
2Ml ± 



m 



.12 



2s^M 2 H± 



(27) 



If 

771^2 is positive and —gs 3> g2i the contributions of the charged Higgs bosons and the W bosons 
interfere constructively in both C and D, and we get an enhancement for large tan/3 or m\ 2 - 
Applying the maximization procedure, one finds that _R 77 ~ 1.7 can be realized for 110 GeV < 
M/j < 150 GeV while still being compatible with all the constraints considered. Regions where this 
may be realized are shown in Fig. [3l 



20.0 



19.8 



19.6 



re 



19.4 



19.2 



19.0 




0.00 



(a) i? 77 with M H = 224 GeV, Ma = 102 GeV and 
mi2 = 50 GeV. Ar w 0.035 in this region and there- 
fore allowed. 




(b) i? 77 with M H = 333 GeV, M A = 111 GeV 
and mi 2 = 100 GeV. Ar w 0.038, which is not 
forbidden at 2 <r. 



Figure 3: analyzed for the V potential (P), for M ft = 130 GeV, M H ± = 100 GeV. 

If m\ 2 is negative, the situation is similar to the Va case as we again get destructive interference 
in both C and D. An enhancement is possible if mf 2 is large, or if tan/3 is large. These regions 
are strongly constrained by the unitarity constraint |a+| < 8ir, which is more restrictive if m\ 2 
is negative. As a+ restricts simultaneous growth of tan/3 and m\ 2 , one can derive a rough lower 
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bound, which is -(250 GeV) 2 < 

771^2 f° r tan /3 ~ 2. Even if such low values were allowed, one 
would only get i? 77 < 0.7 and therefore no enhancement compared to the SM. Lower values of mf 2 
can be attained for smaller values of tan /3, but we would also have to raise M H ± to be compatible 
with the constraint from R^, and therefore this would not result in an overall enhancement either. 

Expanding a+ for large values of tan /3, one sees that the maximal values of tan (3 are smaller for 
non-zero negative values of mf 2 - Performing the numerical maximization, we found no enhance- 
ment, with the partial width staying below the maximal values attained for Va in most of the 
parameter region. 

For the potential Vb we have one free parameter less, as m\ 2 = M\spcp. If M\ ^> M% we 
again have constructive interference in both C and D, and the situation is similar to V with large 
and positive m\ 2 - The maximal enhancement we found is a bit lower, yielding i? 77 ~ 1.6 for 
110 < Mh < 150 GeV 0. Obviously, the parameter space available for such an enhancement is 
smaller for Vb- 

5 Conclusions 

We discussed a possible enhancement of the partial decay width T(h — > 77) in the THDM Type I for 
light charged scalars with M H ± = 100 GeV. Unlike in the Type II model, where M H ± > 300 GeV 
due to constraints from B physics, one cannot exclude light charged scalars in Type I. We maximized 
the ratio of the decay width in the THDM over the decay width in the SM, subject to constraints 
from vacuum stability, tree-level unitarity, Z decays into hadrons, and the Ar parameter. The 
ratio was analyzed for the neutral CP-even Higgs h with 110 GeV < < 150 GeV (which is 
the region where an SM-like light Higgs boson h can be discovered at the LHC in h — > 77), and 
with Mh± = 100 GeV for the charged Higgs boson. The maximal possible enhancement differs 
for the potentials considered due to different interference scenarios. For the IDM the maximal 
enhancement was found to be around +70%, for Va around —20%. The results for V and Vb were 
rather similar, being around 70% and 60%, respectively. This is larger than what was found for 
Type II, where for heavy charged Higgs bosons with M H ± ~ 400 GeV the enhancement was around 



25% [8|. Our results for the enhancement of the partial width in T ype I differ from Ref. 12J, mostly 
because we used the more restrictive unitarity constraints of Ref. [16J and the additional constraint 
R b of Ref. 0. 

The expected accuracy at the LHC for a(pp — > H)B(H — > 77) is around 10 to 15% for an 
integrated luminosity of 100 fb _1 [jj. If the total decay width and the production cross section in 
the THDM Type I are not too different from the SM, the accuracy should be sufficient to distinguish 
between the SM and such THDMs. 

Combining measurements at the photon collider option of the ILC with the e + e~ collider option, 
the partial width of a SM Higgs with = 120 GeV can be determined with 3% accuracy [3]. 
Such measurements, combined with possible direct detection of a charged Higgs may hence be used 
to distinguish between the various possible potentials of the THDM. 



Plots for Vb can be found in Ref. K32 
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